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Abstract 

We compute the coefficients of an infinite family of chiral primary operators in the 
local operator expansion of a circular Wilson loop in M = 4 supersymmetric Yang- 
Mills theory. The computation sums all planar rainbow Feynman graphs. We argue 
that radiative corrections from planar graphs with internal vertices cancel in leading 
orders and we conjecture that they cancel to all orders in perturbation theory. The 
coefficients are non-trivial functions of the 'tHooft coupling and their strong coupling 
limits are in exact agreement with those previously computed using the AdS/CFT 
correspondence. They predict the sub-leading orders in strong coupling and could in 
principle be compared with string theory calculations. 
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1 Introduction and summary of main results 



The idea that a quantized gauge theory could have a dual description as a string theory 
has a long history. Recently one concrete realization of such a duality has emerged. It 
has been conjectured Q that there is an exact mapping between J\f = 4 supersymmetric 
Yang-Mills theory (SYM) with gauge group SU(N) on four dimensional spacetime and 
IIB superstring theory on background AdS§ x S$ with N units of RR flux. 

This mapping is most useful in the low energy, weakly coupled limit of the string 
theory. This coincides with the large N 'tHooft limit [0] of SYM theory, where N 
is taken to infinity holding the combination of Yang-Mills coupling constant and N, 
defined by A = <7 2 N, fixed and then taking the large A limit. This projects onto 
the strong coupling limit of the sum of planar Feynman diagrams. In string theory, 
this coincides with the classical low energy limit where the string theory is accurately 
described by type IIB supergravity on the background space AdS§ x S5. Some explicit 
computations can be done there. The results can then be interpreted in terms of the 
gauge theory using a well-defined prescription M , Q . 

Though it has been used for many computations of the strong coupling limits of 
gauge theory quantities (see refs. |]]-|]] for reviews), it is difficult to obtain a direct 
check of the Maldacena conjecture. The reason for this is the fact that the correspon- 
dence with supergravity computes gauge theory in the large A limit, with corrections 
from tree level string effects being suppressed by powers of 1/y/X and sometimes com- 
putable to the next order. On the other hand, the only other analytical tool which 
can be used systematically in the gauge theory is perturbation theory which is an 
asymptotic expansion in small A. Generally, the only quantities for which these expan- 
sions have an overlapping range of validity is for quantities which are so protected by 
supersymmetry that they do not depend on the coupling constant. 

There is, however, one known example of a quantity which is a non-trivial function 
of the coupling constant and whose large N limit is computable and is thought to be 
known to all orders in perturbation theory in planar diagrams. That quantity is the 



circular Wilson loop. Its expectation value was computed in ref. [1C]. The contribution 
of a subset of all Feynman graphs, the planar rainbow diagrams, were found at each 
order in A and the sum of all orders was taken to obtain the result 

(W [circle]) = -^=h(V\) « as A ^ 00 (1.1) 

where I\{x) is a modified Bessel function. It was also shown explicitly that the leading 
corrections to the sum of rainbow diagrams cancels identically. It was conjectured that 



this cancellation would also occur at higher orders and the result (1.1) was thus the 
exact sum of all planar diagrams. Some support for this conjecture was developed in 
ref. [ 1 1 1 . They also observed that the sum over Feynman diagrams could be obtained 



for all orders in the 1/N expansion and had a beautiful argument that, in the large A 
limit, these higher orders produced the expected higher genus string corrections. 



The large A limit in ( LI ) agrees with the expectation value of the circular Wilson 
loop which was computed using the AdS/CFT correspondence in refs. [jL^l,|i~3|l . If 
( |l . l|) is indeed an exact result, this provides a non-trivial check on the validity of the 
AdS/CFT correspondence. It gives the further interesting possibility of comparing cor- 
rections at sub-leading orders in 1/y/X with string theory computations. Investigations 
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of the relevant string theory technique appears in refs. |]14|1- p9|| but explicit calculations 
of the 1/y/X corrections have not yet been done. 

In this Paper we shall report the computation of a series of expansion coefficients 
which are related to the circular Wilson loop and chiral primary operators in M = 4 
SYM theory. The problem that we pose is the following. When probed from a distance 
much larger than the size of the loop, the Wilson loop operator can be expanded in a 
series of local operators with some coefficients @,||]: 

W[C) = (W[C]) C A R Aa O a {0) (1.2) 
A 

where A (O) is a local operator evaluated at the center of the loop, is the conformal 
dimension of O (x) and R is the radius of the loop. The problem is to compute the 
coefficients Ca in this operator product expansion (OPE). 

We shall concentrate on computing the coefficients for a particular class of chiral 
primary operators (CPO). We will be able to compute the contribution of the sum of 
all planar rainbow graphs to the coefficients Ca in that case. We are also able to show 
that the leading order corrections to this sum, which come from diagrams with internal 
vertices, cancels identically. This leads us to conjecture that the radiative corrections 
cancel to all orders and the sum of planar rainbow graphs gives the exact result. 

We find that the coefficients that we compute are non-trivial functions of the cou- 
pling constant. In the limit of large A they coincide with results of the AdS/CFT 



correspondence [12]. This gives a large array of non-trivial functions of the coupling 
constant which could be compared with string computations of the strong coupling 
limit. There are various reasons why these computations could be simpler than the 
1 /yf~\ corrections to the expectation value of the Wilson loop itself. 

The Wilson loop operator in M = 4 SYM theory that is readily computed using 
the AdS/CFT correspondence and which has the right transformation properties under 
supersymmetry [21], [13] contains the scalar fields inside the path-ordered exponential: 



W[C] = -J- trPexp , 
N IJC 



j dr (iA^xjxn + &{x)6i\x 



(1.3) 



where x^(t) parameterize the contour C and 9i are Cartesian coordinates of a point 
on S 5 : 8 2 = 1. If the size of the contour C is small, the Wilson loop can be expanded 
in local operators as in ( |1.2| ). 

For primary operators, one can choose a basis where 



(0 A {x)0"{y)) = (1.4) 
i x y i 

Then their OPE coefficients can be extracted from the large distance behavior of con- 
nected two-point correlation functions, 



(W(C)0 A {L)) ^ r a a 



Ca TJK7 + --- (1.5) 



(W{C)} A L 2A a 

where L> R and the omitted terms are of higher order in R 2 /L 2 . 
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The coefficient corresponding to the CPO of lowest conformal dimension, which in 
this case is A = 2, is important as it determines the correlator of two Wilson loops 
with large separation 

-d (?)' + ... (1-6) 



(VC[Ci]> (W[C 2 \) '\L 

The coefficients of various CPOs in the expansion of the circular Wilson loop were 
calculated in ref. [|l2| both perturbatively at A ~ and at strong coupling, A ~ 
oo, using the AdS/CFT correspondence. Evaluation of the correlators that define 
coefficients in the strong-coupling regime involves a hybrid of the supergravity and the 
string calculations. 



In AdS/CFT, the Wilson loop operator (1.3) naturally couples to stringy degrees 
of freedom. It creates a classical string world-sheet which is embedded in AdS^ x S§ 



and whose boundary is the contour of the Wilson loop |2l]] , 1 22 ] . 

On the other hand, a local operator O (x) emits one of the supergravity fields 
at point x. When it contributes to a correlator of Oa{x) with the Wilson loop, this 
supergravity mode propagates on the background AdS§ x £5 and is then absorbed by 
a vertex operator which must be integrated over the string world-sheet. 

1.1 Dimension two operators 

Let us begin by considering the CPO with smallest conformal dimension, A = 2. It is 
the symmetric traceless part of a gauge invariant product of scalar fields, 

Qij = tr ( _ I gU^A r 17 \ 

V2X V 6 J K J 

This operator is the lowest weight component of a short multiplet of N = 4 super- 
conformal algebra. Such chiral primary operators have very special properties. The 
super-conformal algebra guarantees that their conformal dimensions do not receive 
radiative corrections, so in this case the conformal dimension is exactly two. Further- 
more, it is known that their two and three-point correlation functions are given by 
the free field values, that is, that they are independent of the coupling constant, g. It 
is known that their four-point functions are non-trivial, so they are not free fields in 
disguise H|-|6|. 

In ([O]), the overall coefficient is chosen to give a canonical normalization of the 
two-point function: 

O ij (x)O kl (y)) = -( 8 ik 6 jl + 5 a 5 jk - - 5 lj 5 kl ) - — l —r . (1.8) 
1 2 V 3 J \x - y| 4 

The small coupling limit of the correlator of O 13 with the Wilson loop is straight- 
forward to obtain. To leading order in perturbation theory: 

(W(C)) -N^ X [" ~6 6 ) U (A ^ 0) - (L9) 

The linear dependence on A is an obvious consequence of the fact that the correlator 
contains two propagators and one power of A is cancelled by the normalization. 
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The AdS dual of the dimension two operator is the negative mass scalar which is a 
linear combination of the trace of the metric and the Ramond-Ramond four-form field. 
Its contribution to the OPE of the circular Wilson loop was calculated in [12]: 

Comparing OPE coefficients at strong and at weak coupling, we see that the scaling 
with A is different. The OPE coefficients are clearly renormalized by radiative correc- 
tions. We shall conjecture that, in the large N limit, this renormalization is entirely 
due to planar rainbow diagrams. We shall also obtain the sum of planar rainbows as 



(W(C)) - N h (yx) 



where I2 and I\ are modified Bessel functions. By construction, this expression reduces 
to (tL9f) at small A. Since 



4 (7a) 

lim — ) f = 1 (1.12) 

for any k, the AdS/CFT prediction ( l.lCj ) is also exactly reproduced at large A. The 
sum of rainbow diagrams thus interpolates between perturbative and strong coupling 
limits of the OPE coefficient. 



1.2 Chiral primary operators 

The O lJ is the first in an infinite sequence of CPOs. The operator of dimension k in 
this sequence is a symmetrized trace of k scalar fields: 

^ = ^^Cj.„ ik tr (1.13) 

where C\ , are totally symmetric traceless tensors normalized as 

cL, k cL, k = 

The AdS duals of CPOs are Kaluza Klein modes of the AdS§ tachyonic scalar on 
S 5 and each CPO is associated with a spherical harmonic: 

Y I (0) = C I il ^ ik e i K..e i K (1.15) 



Here, we are following the notation of refs. [23], [12]. 

The OPE coefficients depend on how the operators are normalized. When compar- 
ing perturbative calculations with the AdS/CFT predictions, we need to use the same 
normalization. For operators in short multiplets of N = 4 supersymmetry, this is easy 
to achieve, since the two point correlation functions of such operators do not receive 
radiative corrections and can be used to fix normalization. The coefficient in ( 1.13j ) is 
chosen to unit normalize the two point function: 

. xU 

I k (x)0{(y)) = w - WJ :. (1.16) 
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The same conventions were used in the supergravity calculations of ref. [|T^] . 

At weak coupling, the OPE coefficient of the circular Wilson loop is proportional 
to A fc / 2 , where k is the dimension of the CPO: 



(W(C)) -N 2 ~kJ L^ Y {9) (A ^ 0) - (L17) 

It turns out that AdS/CFT correspondence predicts a universal scaling of the OPE 
coefficients with A at strong coupling: all of them are proportional to y/\ independently 
of k. This can be easily understood by considering a pair correlator of the Wilson 
loops, which is quadratic in OPE coefficients. The Wilson loop correlator is described 
by an annulus string amplitude and therefore is proportional to the string coupling g s . 
According to the AdS/CFT dictionary, 

9 2 A 



J 4vr 4vriV 

Hence, OPE coefficients must scale as vA- An explicit calculation gives [12]: 



<™> i^VE*^ (1 ,s) 



(W{C)) N L 2k 

Our main result is an expression for correlators of the circular Wilson loop with 
CPOs: 

(w(c) ol) i i k (Vx) R k 

(W(C)) N h L/x) { ' { } 



which we expect is exact in the large N limit. Its expansion in A reproduces ( 1.17 ). The 



strong-coupling limit exactly coincides with the AdS/CFT prediction (using eq. ( 1.1 2| ) ) 



2 Re- summation of rainbow diagrams 

The Euclidean action of M = 4 supersymmetric Yang-Mills theory is 

iV„ fl 

*^ + wr - 

i<j 



J d' Xj Tv -F% + (Z^) 2 - ^'] 2 + 

+t/j T T»D^ - i^T'lV, ^} (2.1) 



where (T^, T % ) are ten dimensional Dirac matrices in the Majorana-Weyl representation. 
We will work in the Feynman gauge where the gauge field propagator has the form 

<<(^)> = A^^^ (2.2) 

Our calculation of the OPE coefficients begins with summing all planar rainbow 
diagrams of the kind shown in fig. They contain k scalar propagators connecting 
the point L to the Wilson loop and propagators of scalars and gauge fields connecting 
points in segments of the loop. 
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Figure 1: A typical diagram that contributes to the correlator of the circular Wilson loop with the 
dimension k CPO. 



If L lies on the axis of symmetry of the circle, the scalar propagators are constants 
equal to 



8vr 2 (L 2 + Ft 2 )' 

If the origin is displaced from the axis of symmetry, the propagators will depend on 
positions of their endpoints on the circle. In any case, we will be interested in the large- 
distance asymptotics, and the propagators can be set to l/8n 2 L 2 , up to corrections of 
higher order in 1/L. 

The problem thus reduces to re-summation of rainbow diagrams for each of the 
segments of the circle. This problem was solved in ref. [27]. In the following we will 
review the salient points involved in finding the solution. We start with the dimension 



two operator (1.7), when there are only two segments. 



2.1 Dimension two operators 

In this case we have: 

( W{C) °") ~ Jr ^tx p? ~ H jf <W*W* - d, 

(2.3) 

where W((p' — <p) denotes the sum of rainbow graphs for a segment of the circle between 
polar angles cp and cp'. To compute this sum, we notice that the sum of the scalar and 
the gluon propagators between any two points on a circle does not depend on the 
positions of these points: 

((M,(,)^ + ^( a; )^|i|)^(U M (y)^ + ^( y )^|y|)J o = p^ Jf^J 

= i^v Wte ' (2 - 4) 
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This observation allows us to replace the field-theory Wick contraction by the matrix- 
model average defined by the partition function 



dM exp 



8^ 
A 



TVtrM 2 



(2.5) 



Upon the replacement of iA^(x)x^ + & % {x)9 % \x\ by M, the sum of rainbow diagrams 
in the segment of the length ip reduces to the matrix-model counterpart of the Wilson 
loop: 



W(cp) 



(Itre^) 



(2.6) 



The matrix model can be viewed as a combinatorial tool which simply counts the 
number of planar graphs [10]. 

The Wilson loop in the matrix model satisfies Schwinger-Dyson identity in the 
large- N limit (the loop equation |28[j): 



A 



16vr 2 



dtpW{<p)W(#-y). 



The solution |1CJ , |2{J , [ 27J of the loop equation is 



W(<p) 



47T 



h 



2tt 



(2.7) 



{21 



The integral we need to compute in order to calculate the correlation function ( [2.3] ) 
is the right-hand-side of the loop equation at 6 = 2tt. Using properties of the modified 
Bessel functions, 



kl k {z) 



{h-i(z) +h+i{z)) , 
(h-i(z) - I k +i(z)) , 



we get: 



dyW{ip)W{-d-y) 



32vr 2 ( y/X-d y 
h 



2tt 



Setting 9 = 2tt and substituting into (|2.3[) , we obtain: 



W{C) O i 



^2V2/ 2 (Va' 



*Q3 _ 

6 



L 4 



Dividing by the vacuum expectation value of the Wilson loop, 

(W(C)) = W(2tt) = ^=h(VX) , 



(2.9) 
(2.10) 

(2.11) 

(2.12) 



we arrive at the result (1.11) which we quoted earlier. 
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2.2 Chiral primary operators 

The correlator of the Wilson loop with the CPO of dimension k contains an integral 
over k — 1 endpoints of the scalar propagators (one integration yields an overall factor 
of 2vr): 



W{C)0{) = - /Q _ 2 r2 , fc C^jfr ■ ■ R h 



1 (87T 2 ) k / 2 

N ^k\k/2 (8vr 2 L 2 )* 

p27T C i Pk — 2 

x2tt / dept... difk^ W(y>fc_i)W(v fc _2 - tPk-i) ■ ■ ■ W(2ir 

Jo Jo 

It is useful to introduce 

rip f L Pk — 2 

F k (cp)= di Pl ... dtp k ^iW((p k - 1 )W((pk-2-<Pk-i)-''W(ip-(pi). (2.14) 

Jo Jo 

The correlator is expressed in terms of F^^tt) as 

To find the functions Fk((p), we again use the loop equation. Differentiating Fk(^p) 
and using (0), we get the recurrence relations: 

Ftfr) = F k -i(<p) + J^2 F k+i(<P), (2-16) 

which are supplemented by initial conditions 

F 1 (tp) = W(<p), F ((p)=0. (2.17) 

These unambiguously determine all F k . A systematic way to solve these recurrence 
relations is to introduce a generating function and then use a Laplace transform to 
convert differential equations into algebraic equations. The result is 

It is straightforward to check that this expression solves the recurrence relations with 
the help of Q. 

Substituting ( gig) into ( |2~T^ ), we obtain 

W(C) Oi) = 1 2 k ' 2 VkI k (VX) Y\6). (2.19) 



N V / L 2fc 

Normalizing by the Wilson loop expectation value (|2.12|) , we get ( |1.19|) . 



3 Radiative Corrections 

The leading radiative corrections come form the Feynman diagrams which are shown 
in figgj. 

*The diagram similar to e, but with scalar lines replaced by gluon propagators does not contribute because 
of R-charge conservation. 
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Figure 2: Leading radiative corrections to the correlator of the Wilson loop operator with the 
k = 2 CPO. These are not taken into account by the sum over planar rainbow graphs which were 
computed in Sec. 3. 



Each of these diagrams is separately divergent and regularization is required to 
define them properly. We use a regularization by dimensional reduction which was 
previously used in ref. pp| . The essential observation is that Af = 4 SYM is obtained 
by dimensional reduction of ten dimensional SYM. This dimensional reduction retains 
sixteen supersymmetries in any dimension. Thus, a supersymmetric dimensional regu- 
larization of N = 4 SYM theory is obtained by dimensionally reducing ten dimensions 
to 4 — e dimensions. 

In this dimensional regularization, the diagram in fig.^a is of higher order than the 
relevant leading power, R 2 /L 4 , and therefore does not contribute to C^. 

In the limit L S> R, using dimensional regularization, the leading, R 2 /L 4 , contri- 
butions of the remaining diagrams in fig.|| can be seen to be identical to results of 
computing the diagrams which are displayed in fig.||. This sum of diagrams is known 




Figure 3: The contribution to the leading term in R 2 /L 4 of the diagrams in fig.|2|b-e are given 
by these 2-loop diagrams. This combination of diagrams is known to vanish due to the non- 
renormalization theorem of the 2-point function of the CPO. 

to vanish when the dimension is exactly four, due to the non-renormalization the- 
orem for the two-point function of the CPO. This non-renormalization results from 
super-conformal invar iance. 
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Similar arguments apply to the higher CPO's for which similar non-renormalization 
theorems can be applied |3Cj] . 

This is by no means a proof that all radiative corrections vanish. But the excellent 
agreement with strong coupling AdS / CFT results gives optimism that it is indeed the 
case. 



4 Remarks 



Our main results ( |1.11| ) and ( 1.1 9| ) are valid when the distance from operator insertion 



to the loop is much larger than the loop's radius. However, it is not hard to restore 
the dependence on the radius and on the orientation of the circle. Consider first the 
case when the operator is inserted at the symmetry axis of the circle. As follows from 
the discussion in the beginning of Sec. g we can find the correlators of the Wilson loop 
with CPOs at any R and L, not only at large distances, simply by replacing I? with 
L 2 + R 2 . As expected, this is perfectly consistent with the supergravity prediction, 
which actually is known in complete generality. If we denote the displacement of the 
operator insertion from the axis of symmetry by r and, as before, L is the distance 
from the plane of the circle, the correlation function is obtained from the large-distance 
asymptotics by replacing 

^ (41) 

L 2k [(L 2 + r 2 -R 2 ) 2 + AL 2 R 2 ] k/2 

At r = 0, the denominator on the right hand side indeed coincides with {I? + R 2 ) k . 

We were not able to compute the sum of rainbow diagrams for but the validity 

of the prescription ([O]) can be proved by an indirect argument. The point is that, 
once we know the correlators at r = 0, the dependence on r is unambiguously fixed 
by conformal invariance, because r can be always set to zero by a special conformal 
transformation, which maps a circle onto a circle. This transformation is not anomalous 
and therefore does not affect the dependence of the correlator on A. 

Another remark concerns the dependence of the OPE coefficients on A:. It was 
argued on general grounds that coupling of a Wilson loop to states of very large spin 



should be factorially suppressed [31]. In a theory with an exponential density of states, 
unsuppressed coupling to such states would lead to a catastrophe in the pair correlator 
of Wilson loops similar to the Hagedorn transition. The CPOs we consider in this paper 
carry the spin of £0(6) R-symmetry which for 0\ is equal to k. At weak coupling, 
the OPE coefficients are indeed suppressed at large k, as follows from ( 1.17| ), but the 



supergravity result ( l.lEj ) seems to suggest that this suppression disappears at strong 



coupling. Careful inspection of the exact OPE coefficients ( 1.1S| ) shows that limits 



A — » oo and k — > oo do not commute and coupling of the Wilson loop to operators 
of very high spin is always suppressed, but the suppression begins with operators of 
parametrically large spin k ~ A. At k S> A: 



fc! h (yx) • ( J 



n 



5 Discussion 



We have calculated the OPE coefficients of the circular Wilson loops by re-summation 
of the planar Feynman graphs without internal vertices. It is likely that other diagrams 
cancel to all orders of perturbation theory. We have successfully checked this conjecture 
up to two loops. Complete agreement of an infinite set of OPE coefficients with the 
supergravity predictions at strong coupling strongly suggests that our results are indeed 



exact in the 't Hooft limit. It would be interesting to see if the arguments of Ref. [11| 
based on diagram-by-diagram conformal transformations can be invoked to prove that 
only rainbow diagrams contribute to the OPE coefficients of the circular Wilson loop 
with CPOs. 

Our result for the OPE coefficients, along with the known exact expression for the 
vacuum expectation value of the circular loop, can be regarded as a prediction for the 
string theory in AdS§ x S 5 . The usual a' expansion of the world-sheet sigma-model 
then coincides with the expansion in l/\/A: 



^ C) °V _ 1 a*/*- W (l - + fc4 - 4fc2 + 3 + ( J_V\ *L YHO) 

(W(C)} ~ N VfcA l v 1 2^A 8A +U [ X 3/2)) L 2k 

(5-1) 

The calculation of the stringy correction to the expectation of the Wilson loop is a 
hard problem analogous to instanton calculations in field theory (see Jl^j, for details). 
As usual, such problems require delicate treatment of various normalization factors 
associated with zero modes and with regularization of fluctuation determinants. How- 



ever, in the ratio of the expectation values (5.1), these normalization factors cancel. 
For this reason, a calculation of stringy corrections to the OPE coefficients seems less 
complicated and perhaps can be accomplished without tremendous effort. 

It would also be interesting to consider similar correlators of Wilson loops with 
different contours [[O] or with other operators, where some preliminary results are 



contained in ref. [33], [34] 



Acknowledgments 

We are grateful to Jan Plefka, Matthias Staudacher and Arkady Tseytlin for useful 
comments. This work was supported by NSERC of Canada and NATO Collaborative 
Linkage Grant SA(PST.CLG. 977361)5941. The work of K.Z. was also supported in 
part by the Pacific Institute for the Mathematical Sciences and in part by RFBR grant 
01-01-00549 and RFBR grant 00-15-96557 for the promotion of scientific schools. 



References 

[1] J. Maldacena, "The large N limit of super-conformal field theories and supergrav- 
ity," Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 (1998)] 
||hep-th/9711200|| 



[2] G. 't Hooft, "A Planar Diagram Theory For Strong Interactions," Nucl. Phys. B 
72, 461 (1974). 

[3] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, "Gauge theory correlators from 



non-critical string theory," Phys. Lett. B 428, 105 (1998) jhep-th/ 9802109 1. 



12 



E. Witten, "Anti-de Sitter space and holography," Adv. Theor. Math. Phys. 2, 
253 (1998) |hep-th/9802150| . 

J. L. Petersen, "Introduction to the Maldacena conjecture on AdS/CFT," Int. J. 
Mod. Phys. A 14, 3597 (1999) jhep-th/9902T3lH . 

P. Di Vecchia, "Large N gauge theories and AdS/CFT correspondence," hep 



th/9908148 



E. T. Akhmedov, "Introduction to the AdS/CFT correspondence," hep- 



th/9911095 



O. Aharony, S. S. Gubser, J. Maldacena, H. Ooguri and Y. Oz, "Large N field 
theories, string theory and gravity," Phys. Rept. 323, 183 (2000) jhep-th/9905111|1 . 
D. Z. Freedman and P. Henry-Labordere, "Field theory insight from the AdS/CFT 



correspondence," hep-th/001108(: . 

J. K. Erickson, G. W. Semenoff and K. Zarembo, "Wilson loops in N 



super- 



symmetric Yang-Mills theory," Nucl. Phys. B 582, 155 (2000) [|hep-th/0003055 l. 
N. Drukker and D. J. Gross, "An exact prediction of N = 4 SUSYM theory for 



string theory," |hep-th/0010274 . 

D. Berenstein, R. Corrado, W. Fischler and J. Maldacena, "The operator product 
expansion for Wilson loops and surfaces in the large N limit," Phys. Rev. D 59, 
105023 (1999) [|hep-th/9809T88|| . 

N. Drukker, D. J. Gross and H. Ooguri, "Wilson loops and minimal surfaces," 
Phys. Rev. D 60, 125006 (1999) jhep-th/9904191[| . 

J. Greensite and P. Olesen, "Worldsheet fluctuations and the heavy quark poten- 
tial in the AdS/CFT approach," JHEP 9904, 001 (1999) jiep-th/9901057H 
S. Forste, D. Ghoshal and S. Theisen, "Wilson loop via AdS/CFT duality, 



hep- 



th/0003068 



S. Forste, D. Ghoshal and S. Theisen, "Wilson loop in N = 4 super Yang-Mills 
theory," Fortsch. Phys. 48, 111 (2000). 

S. Forste, D. Ghoshal and S. Theisen, "Stringy corrections to the Wilson loop in 

N = 4 super Yang-Mills theory," JHEP 9908, 013 (1999) [|iep-th/9903042l . 

Y. Kinar, E. Schreiber, J. Sonnenschein and N. Weiss, "Quantum fluctuations of 



Wilson loops from string models," Nucl. Phys. B 583, 76 (2000) [|hep-th/99111~23l 
N. Drukker, D. J. Gross and A. Tseytlin, "Green-Schwarz string in AdS(5) x S(5): 
Semiclassical partition function," JHEP 0004, 021 (2000) [|hep-th/0001204]| . 
M. A. Shifman, "Wilson Loop In Vacuum Fields," Nucl. Phys. B 173, 13 (1980). 
J. Maldacena, "Wilson loops in large N field theories," Phys. Rev. Lett. 80, 4859 
(1998) lhep-th/9803002l . 

S. Rey and J. Yee, "Macroscopic strings as heavy quarks in large N gauge theory 



and anti-de Sitter supergravity," [hep-th/9803001 . 

S. Lee, S. Minwalla, M. Rangamani and N. Seiberg, "Three-point functions of 
chiral operators in D = 4, N = 4 SYM at large N," Adv. Theor. Math. Phys. 2, 



697 (1998) |hep-th/9806074 |. 

G. Arutyunov, B. Eden, A. C. Petkou and E. Sokatchev, "Exceptional non- 
renormalization properties and OPE analysis of chiral four-point functions in N 



= 4 SYM(4)," |hep-th/010323C . 

L. Hoffmann, A. C. Petkou and W. Ruhl, "Aspects of the conformal operator 



product expansion in AdS/CFT correspondence," [hep-th/0002T54] . 

M. Bianchi, S. Kovacs, G. Rossi and Y. S. Stanev, "On the logarithmic behavior 

in N = 4 SYM theory," JHEP 9908, 020 (1999) [|hep-th/9906188|l . 



13 



[27] K. Zarembo, "String breaking from ladder diagrams in SYM theory," JHEP 0103, 

042 (2001) ||hep-th/0103058| . 
[28] A. A. Migdal, "Loop Equations And 1/N Expansion," Phys. Rept. 102, 199 

(1983). 

[29] G. Akemann and P. H. Damgaard, "Wilson loops in N = 4 supersymmetric Yang- 
Mills theory from random matrix theory," hcp-th/0101225 . 

[30] E. D'Hoker, D. Z. Freedman and W. Skiba, "Field theory tests for correlators in 
the AdS/CFT correspondence," Phys. Rev. D 59, 045008 (1999) [|hep-th/980709l| . 

[31] K. Zarembo, "Wilson loop correlator in the AdS/CFT correspondence," Phys. 
Lett. B 459, 527 (1999) flhep-th/9904149] . 

[32] J. K. Erickson, G. W. Semenoff, R. J. Szabo and K. Zarembo, "Static potential 
in N = 4 supersymmetric Yang-Mills theory," Phys. Rev. D 61, 105006 (2000) 
| hcp-th/991 1088|| . 

[33] U. H. Danielsson, E. Keski-Vakkuri and M. Kruczenski, "Vacua, propagators, and 

holographic probes in AdS/CFT," JHEP 9901, 002 (1999) [|hep-th / 981 2007 ]. 
[34] J. Erickson, G. W. Semenoff and K. Zarembo, "BPS vs. non-BPS Wilson loops in 
N = 4 supersymmetric Yang-Mills theory," Phys. Lett. B 466, 239 (1999) flh"ep 



th/9906211 ]. 



14 



